Explicit equilibria in a kinetic model of gambling 
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We introduce and discuss a nonlinear kinetic equation of Boltzmann type which describes the evolution of 
wealth in a pure gambling process, where the entire sum of wealths of two agents is up for gambling, and ran- 
domly shared between the agents. For this equation the analytical form of the steady states is found for various 
realizations of the random fraction of the sum which is shared to the agents. Among others, Gibbs distribution 
appears as steady state in case of a uniformly distributed random fraction, while Gamma distribution appears 
for a random fraction which is Beta distributed. The case in which the gambling game is only conservative-in- 
the-mean is shown to lead to an explicit heavy tailed distribution. 



I. INTRODUCTION 

Various concepts and techniques of statistical mechanics 
have been fruitfully applied for years to a wide variety of com- 
plex extended systems, physical and otherwise, in an effort to 
understand the emergent properties appearing in them. Eco- 
nomics is, by far, one of the complex extended systems to 
which methods borrowed from statistical mechanics for parti- 
cle systems have been applied ||7lll[TTl[T3l[Tl|23l|24i|29ll. In 
most of the models introduced so far, the trading mechanism 
leaves the total mean wealth unchanged. Then, a substantial 
difference on the final behavior of the model (presence or not 
of tailed steady states) can be observed depending on the fact 
that binary trades are pointwise conservative, or conservative 
in the mean lfT5ll2^ . The asymptotic distribution of wealth, 
however, depends completely on the microscopic structure of 
binary trades. Other kinetic models have been recently pro- 
posed, which, while maintaining the kinetic description, in- 
troduce more sophisticated rules for trading. For example, a 
description of the behavior of a stock pric e has been developed 
by Cordier, Pareschi and Piatecki in [|l2|l ■ Further, there have 
been efforts to include non-microscopic effects, like global 
taxation (and subsequent redistribution), in recent works of 
Guala [22], Pianegonda, Iglesias, Abramson and Vega ll28ll . 
Garibaldi, Scalas and Viarengo ll9ll and Bisi, Spiga and the 
present author B IsTll . 

Despite the high number of studies devoted to the sub- 
ject, well documented by various recent review papers ^ |2ll 
^tIHSHH, analytical solutions or explicit steady states for 
wealth distribution densities are rarely present in the litera- 
ture. The few exceptions are the self-similar solution with 
Pareto tails found by Slanina fl^ for a kinetic model of a 
non-conservative (decreasing in the mean) economy, and the 
steady state solution (Gibbs distribution) of a linear kinetic 
equation modelling taxation and uniform redistribution lIsTll . 

The goal of this paper is to show that in a pure gambling 
trade market explicit equilibrium solutions can be obtained by 
suitably choosing the random fraction of money which gov- 



erns the microscopic interaction. 

Wealth exchange processes in a discrete market are char- 
acterized by binary trades. A fixed number of agents in a 
system are allowed to interact (trade) stochastically and thus 
wealth is exchanged between them. The basic steps of such a 
wealth exchange model can be described as follows 

Wi(f + 1) = VV;(f)+Aw, Wj{t+\)^Wj{t)-Aw, (1.1) 

where w,(f) and Wj{t) are wealths of /-th and y'-th agents at 
time t and Wi{t + 1) and Wj{t +1) are that at the next time 
step f + 1 . The amount Aw (to be won or to be lost by an 
agent) is determined by the nature of interaction. If the agents 
are allowed to interact for a long enough time, a steady state 
equilibrium distribution for individual wealth is achieved. The 
equilibrium distribution does not depend on the initial config- 
uration (initial distribution of wealth among the agents). 

In a pure gambling process [ 14], the entire sum of wealths 
of two agents is up for gambling. Some random fraction of 
this sum is shared by one agent and the rest goes to the other 
The randomness is introduced into the model through a pa- 
rameter e which is a random number drawn from a probability 
distribution in [0, 1]. In general it is assumed that e is indepen- 
dent of a pair of agents, so that a pair of agents is not likely 
to share the same fraction of aggregate wealth the same way 
when they interact repeatedly. The interaction rule can be seen 
through 



wi{t + \) = e[wi{t) + wj{t)], 
wj{t + l) = {l-e)[wi{t) + wj{t)], 



(1.2) 
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where the pair of agents (indicated by i and j) are chosen 
randomly. The amount of wealth that is exchanged is now 
Aw = e[w,(f) +Wj(t)] — w,(f). Numerical experiments fl^ 
show that, if e is a random number drawn from a uniform dis- 
tribution in [0,1], the individual wealth distribution at equi- 
librium emerges out to be Gibbs distribution like exponential. 
It is strongly believed, and somewhat argued from numerical 
experiments, that a slight variant of the previous model pro- 
duces at equilibrium Gamma-like distributions [10], but, as 
remarked by the authors, the form of the exact solution is still 
an open question. The main feature of trade ( 11.21 ) is such the 
amount of money which is restituted in a single trade coin- 
cides with the amount of money the two agents use for gam- 
bling {pointwise conservative trade). In this gambling market 
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there is no rule of the (conserved) total amount of money ini- 
tially in the hands of the agents. In other words, agents do 
not take advantage from the amount of money available in the 
market. 

The idea of using this money as a reservoir for trades can 
be easily implemented allowing agents to trade with random 
profit 



wi{t + l) = ex[wi{t) + wj{t)l 
Wj{t + \) = e2[wi{t) + Wj{t)], 



(1.3) 



where the parameters (ei , £2) are now random numbers drawn 
from a joint probability distribution such that 



(ei+e2) = 1, 



(1.4) 



where (•) denotes as usual the mathematical expectation. 
Within this picture, Wj(f + 1) + Wj{t + 1) can be strictly less 
than Wi{t) + Wj{t), and in this case the lost money is achieved 
by the market, or the reverse situation is verified, and the addi- 
tional money is taken from the market. Condition ( 11.41 ). how- 
ever, guarantees that in the mean the wealth present in the mar- 
ket is left unchanged. Note that if one assumes that > 7 > 0, 
i— 1,2 then agents are prevented from loosing all their money 
in a single trade. Note also that, choosing (ci , £2) = (e, 1 — e) 
the collision rule ( ll.3l l becomes (11.2b . 

In order to produce a fair game, it will be assumed that 
the random numbers £1 and £2 are identically distributed. For 
trades of type ( ll.2l i. this assumption simply forces the random 
fraction £ to be symmetric with respect to the value 1 /2. 

At a continuous level, the evolution of the wealth distribu- 
tion can be described by means of a nonlinear kinetic model of 
Boltzmann type, in which the wealth distribution of the mar- 
ket is driven by collisions of type ( 11.21 ) |6]. Since the total 
mean wealth is maintained constant in time, the wealth distri- 
bution is known to converge exponentially towards a steady 
profile ifTsl l26[l , which depends on the details of the trade 
mechanism through £,. It is well-known that the continuous 
description takes advantage from the the possibility to make 
use of the Laplace transform version of the kinetic equation 

a. 

Owing to this representation, we will prove that various 
cases are explicitly solvable, and lead to an analytic expres- 
sion of the steady profile. In particular, choosing £, to be sym- 
metric Beta distributions one obtains as equilibrium a Gamma 
distribution with a parameter which depends on the parameter 
of the Beta distribution. A special case emerges here, since 
Gibbs distribution emerges for £,■ uniformly distributed. Also, 
the case of the winner takes all game can be studied in details 
as limit of Beta £, with parameter tending to zero. In all these 
cases, however, the equilibria possess moments of any order 

Interestingly enough, the case of the conservative in the 
mean gambling trade ( |1.3t can be treated likewise, by choos- 
ing the random variables £,, / = 1 , 2 to be inverse Beta. In this 
case, however, the steady state distribution is shown to be a 
generalized Gamma distribution with fat tails, which contains 
as a particular case the distribution found by Slanina |29] in 
a different context. This founding clarifies through a simple 
example the role of the social use of the money present in the 



market. The possibility to access to the money available in the 
common reservoir allows the formation of the rich class. 

The paper is organized as follows. In the next section we 
introduce the continuous model which is described by a non- 
linear kinetic equation of Boltzmann type, and its main fea- 
tures are discussed in some detail. The analytical solutions 
in the case of the pure gambling trade ( 11.21 ) are described in 
Section 3. Section 4 deals with the gambhng rule (11.3b . 



II. A CONTINUOUS KINETIC MODEL 

Given a fixed number of agents in a system, which are 
allowed to trade, the interaction rules (11. 3b describe a stochas- 
tic process of the vector variable (vi (f), . . . , V7v(/)) in discrete 
time t. Processes of this type are thoroughly studied e. g. 
in the context of kinetic theory of ideal gases. Indeed, if the 
variables v, are interpreted as energies corresponding to /-th 
particle, one can map the process to the mean-field limit of the 
Maxwell model of elastic particles lUll^. The full information 
about the process in time t is contained in the A^-particle joint 
probability distribution PA?(f , vi , V2, . . . , v^f). However, one can 
write a kinetic equation for one-marginal distribution function 

Pl{t,v) = j PN{t,V,V2,...,VN)dV2---dVN, 

involving only one- and two-particle distribution functions llH 
i. 

Pi (f + 1 , y)-Pi (f , y) = [| P2(f , V,-, y,-) (5(y - £1 (y,- + y,)) 
+ 5(y - £2) (y,- + Vj))yv, dvj - 2Pi (f , y)] ^ , 

which may be continued to give eventually an infinite hierar- 
chy of equations of BBGKY type Q. The standard approxi- 
mation, which neglects the correlations between the wealth of 
the agents induced by the trade gives the factorization 

P2it,Vi,Vj)^Pl{t,Vi)Pl{t,Vj), 

which implies a closure of the hierarchy at the lowest level. In 
fact, this approximation becomes exact for 00. Therefore, 
in thermodynamic limit the one-particle distribution function 
bears all information. Rescaling the time as x = 2t /N in the 
thermodynamic limit N ^ °°, one obtains for the one-particle 
distribution function f{v,t) = Pi(v,t) the Boltzmann-like ki- 
netic equation 



dt 



\{j f{t,^i)f{t,vj)(5{v-e,(vi + v,)) 

+ 5{v-e2){vi + Vj))^dVidv^~f{t,v), 



(2.1) 



which describes the process ( 11.31 ) in the limit N ^ °°. 

Owing to equations of type (12.1b . the study of the time- 
evolution of the wealth distribution among individuals in a 
simple economy, together with a reasonable explanation of 
the formation of tails in this distribution has been recently 
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achieved in 1*2^ (see also lfT5l[l6ll ). The Bohzmann-Hke equa- 
tion ( 12. It can be fruitfully written in weak form. It corre- 
sponds to say that the solution /(v,f) satisfies, for all smooth 
functions ^(v) 



(2.2) 



- (j) (v) - <j) {w))f{t,v)f{t,w) dvdw'j , 

where the post-trade wealths (v* , w* ) obey to the rule ( 11.31 1 

v*=£i(v + w), w*=£2(v + w). (2.3) 

Note that (|2.2| i implies that f{v,t) remains a probabihty den- 
sity if it so initially 



fit,v)dv^ [ Mv)dv=l. 



(2.4) 



Moreover, on the basis of (ll.4l i. the choice (j){v) ~ v shows 
that also the total mean wealth is preserved in time 

m{t)^ [ vf{t,v)dv^ [ vfo{v)dv^m{0). (2.5) 

Consequently, without loss of generality, in what follows, we 
assign to the initial density a unit mean 



/ v'/o(v)c/v=l. 



(2.6) 



Setting 0(v') = exp{-^v} in (|Z2] |. H], one gets the Boltz- 
mann equation for the Laplace transform / of /, where 



Direct computations |26] show that f{t,£,) satisfies the equa- 
tion 

^^^+fit,^) = l{fit,^e,)'+f{t,^e2?). (2.7) 

As extensively discussed in ITsl l26ll . explicitly computable 
conditions on the gambling variables i — 1,2, guarantee 
that the distribution f{t,v)dv converges weakly to a universal 
probability distribution whose Laplace transform is the unique 
solution of 

M^) = \{M^eif+M^e2)^) (2.8) 



£(<^)l^=o = -i. 



with 



This fact follows from Thm. 3.3 in 112611 (see also Thm. 2 in 
which links both the convergence and the boundedness 
of moments of the equilibrium solution to the sign of the key 
function G{s), defined as 



In particular the convergence result is valid without additional 
assumption on /o provided that G '(1) < 0. 

By further requiring that the random variables Ci and £2 are 
distributed with the same law, the equilibrium solution /„ is 
found to be the unique solution to 



(2.9) 



where the random variable £ is distributed according to the 
common law of £1 and £2. Equation ( 12. 9t can be better under- 
stood by saying that, if Z is a random variable with law f^, 
the law of Z, defined in ( 12.91 ), is a distributional fixed point of 
the equation 



Z=''£(Zi+Z2), 



(2.10) 



where —'^ means identity in distribution and one assumes that 
the random variables Zi,Z2 and Z have the same probabil- 
ity law, while the variables Zi,Z2 and £ are assumed to be 
stochastically independent. Equations of ty pe (12.101 ) are well 
known and extensively studied, see e.g. 



ty pe (jz, 

mm- 



III. THE PURE GAMBLING TRADE MODEL 

Let us start from the pure gambhng case in which (£1 , £2) = 
(£,! — £). Having in mind that numerical experiments are usu- 
ally done with a random number drawn from a uniform dis- 
tribution in [0, 1], which leads the discrete market to a Gibbs 
distribution at equilibrium [21], we fix the random number £ 
to be a symmetric Beta random variable of parameters {a, a), 
with (fl > 0). We recall that a random variable is Beta dis- 
tributed with parameters {a,b) if its density is 

= ^J^x-Hi-^r^ xe (0,1), 
r{a)r{b) 

see, for instance, BlSll . 

The case a = I, where £ is a random number uniformly 
distributed on (0, 1), confirms the numerical outcome. In this 
case, in fact, ( 12.91 ) becomes 



/ mx)dx. 

Jo 



It is easy to see that 

is a solution of i2.9i . such that /cI(^)|^=o ~ Since (1 + 
^ ) ^ ' is the Laplace transform of the Gibbs distribution of unit 



(v>0), 



G(5):=(£j + £|)-1 



Gibbs distribution results as analytical steady solution to the 
pure gambling trade market, in case £ is a uniform random 
number in (0, 1). 
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To treat the more general case it suffices to recall that, if 
ai,a2,^ > 0, then formula 3 . 1 97 .4 in [20] gives 



r(ai + a2) 



r(ai)r(a2)Jo (l + ^x)«i+«2 
= (l + <^)-"'. 



dx 



(3.11) 



Using identity (13.1 11 1 with ai = a2 = a, one obtains that for a 
general a > the function 



(3.12) 



solves equation ( 12.91 ). Hence, the equilibrium solution results 
in a Gamma distribution of unit mean, with shape parameter 
a and scale parameter 1/a 



/»(v) = 



r(a) 



(v>0). 



(3.13) 



To verify that the results of [26] and HI], which guarantee 
the convergence to the steady state together with its mo- 
ment boundedness properties hold, one needs to show that 
G'(l) <0 holds true for 

G(^):=2(e-^-)-l, 

where £ is a symmetric Beta random variable of parameters 
{a, a), with (a > 0). In this case, however, 

G'(l) (eloge) <0 

is a simple consequence of the fact that < e < 1 with prob- 
ability one. 

The uniform distribution (a — I) appears like a natural sep- 
aration between two different behaviors of the equilibrium so- 
lutions. In case a < 1, Gamma distributions ( 13.131 1 are mono- 
tonically decreasing, starting from/oc(0) = +°°. In the oppo- 
site case a > 1, there is appearance of a peak around the unit 
mean value. The average wealth is unchanged for every a but, 
when a > 1 , the number of agents with a wealth closer to the 
average value increases or, in other words, the wealth distribu- 
tion becomes more fair for larger a. Analogously, for a < 1 the 
distribution gives more weight to agents with a wealth close 
to zero. 

As a consequence, measures of the inequality of the wealth 
distribution, such as the Gini coefficient, increase for decreas- 
ing a and tend to zero for a Some insight can be gained 
by the simple computation of the variance of /„o. It holds 



VariU) = / v\U{v)d 



v-l^i. 
a 



Hence the variance (the spreading) decreases as a increases. 

There are two interesting limiting cases. The first one is 
obtained by letting a In this case £ converges in dis- 

tribution to the constant value 1/2 and one can immediately 
see that the steady state is a degenerate distribution concen- 
trated on the value 1 of the mean wealth. This corresponds to 
a perfectly fair distribution, in which all agents end up with 



the same wealth. The other limit case is the winner takes all 
game, which can be obtained letting a — > 0. In this case the 
steady state is concentrating on 0, while its variance is blow- 
ing up. This corresponds to the discrete situation in which 
a finite number of agents end up with no money, except one 
which takes all. In the continuous case, the value a — can 
not be assumed directly, since the exchange of the limits a^Q 
and f oo is not allowed by the lack of regularity of the equi- 
librium solution for a ~Q. 



IV. ANALYTIC EQUILIBRIA WITH HEAVY TAILS 

Let us now consider the gambUng rule il.3i . In agreement 
with the previous section, we assume that 



£; — — , i — 1,2, 
' 40,' ' ' 



(4.14) 



where, for a given a > 1, 0,- /= 1,2 is a Beta random variable 
of parameters (a + 1 /2,a— 1 /2). It is immediate to reckon 
that the random variables £,■ are such that <£,■>= 1/2 /= 1,2. 
Consequently (12.5b holds and the model is conservative in the 
mean. We can invoke again Thm. 3.3 in to prove that 
there is a unique equilibrium solution such that its Laplace 
transform satisfies 



U^) = U^/Mflia+l/2,c,-l/2{x)dx, 



(4.15) 



with /oo(^)|^=o = ^1. In this case to apply the results of |2^ 
and |1], which guarantee the convergence to the steady state 
together with its moment boundedness properties, one needs 
to show that G ' ( 1 ) < holds true for 



G(.):=(e[ + e^>-l 
/■' 1 



(4.16) 



2'-^T(2fl)r(fl-5 + i 

r(2a-.?)r(fl + i) 



1. 



The proof of this condition is not direct. For the sake of 
brevity, we postpone the computations to the Appendix. 

Now, we shall prove that the solution of (14.15b is obtained 
by taking the Laplace transform of the so called Inverse- 
Gamma distribution [30] of shape parameter a and scale pa- 
rameter a — I, that is 



(a-l)«e 



_{a-l) 



r{a) 



(4.17) 



which is peaked around the mean value 1 and has heavy tails, 
in that it decays at infinity like An analytical solu- 

tion of type ( 14.17b . with a polynomial decay corresponding to 
fl = 3/2 has been discovered by Slanina [29] as self-similar 
solution of a kinetic model of a non-conservative (decreasing 
in the mean) economy. Motivated by the analogy with a dis- 
sipative Maxwell gas, in 1.29.] a decreasing in the mean wealth 
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model where 



V = pv + qw, w = qv + pw; 
p>q>0, + 



(4.18) 



has been discussed. Note that, within condition ( 14.181 1 on the 
mixing parameters p and q, 

V* +w* = (1 ~2y/pq) (v + w) < v + w, 

which implies that the mean value m(t) at time f decays expo- 
nentially to zero at the rate 1^/pq. The standard way to look 
for self-similarity is to scale the solution. More precisely, de- 
fine the rescaled solution g by 



Kf,v) ==m(f)/(f,m(f)v) 



which imphes that / vg(f , v) dv = I for all f > 0. 

In terms of the Laplace transform g of g, it is found that the 
equation satisfied by g reads I I 6.1 



^ + ^{p + q-l)^=g{p^)g(q^)-g{^] 
Steady solutions to equation (14.20b satisfy 

^ip+q-i)^=m)m)-m- 

Direct computations then show that the function 



(4.20) 



(4.21) 



(4.22) 



solves ( 14.211 1 for all values of p and q satisfying the constraint 
y/p + y/q = 1 . Note that ( I4.22l i is the explicit Laplace trans- 
form of 



4v) = 



(1/2)3/2^-27 

r(3/2) vV2 



(4.23) 



Let us set p = q = 1/4 in ( 14.211 ). Then the steady solution 
( |4.22| l satisfies 



2 



(4.24) 



Following lH, Sect. 6, KM can be equivalently written in 
integral form as 



dp, 



or, setting p'/^ = ^ 



— I xdx, 
4x , 



(4.25) 



(4.26) 



which is nothing but ( 14.151 ) with a = 3/2. Consequently the 
distribution ( |4.23l l solves ( 14.151 ) with a — 3/2. This argument 
establishes a connection between the present problem and the 



non-conservative one introduced by Slanina ||29|] . which leads 
to explicit computations. 

In order to prove that the Laplace transform of ( 14.171 ) is 
the solution of ( 14.151 1. since a direct prove seems not straight- 
forward like in the previous case, we recast the problem in a 
more probabilistic way. First of all, let us note that an Inverse- 
Gamma random variable Y of parameter — 1) can be ob- 
tained by taking Y = 1/X where X has Gamma distribution of 
parameter (a, l/(fl — 1)). Recall that X is a random variable 
Gamma {a, I /{a — 1)) if its density is 



{a-iyv''-^e 



l„(a-l)v 



(4. 19) Recall also that its Laplace transform reads 



(l+^/(«-l))-«. 



(4.27) 



As discussed in Section|lll equation ( 14.15b can be rewritten in 
equivalent way as 



1 

40 



[Y1+Y2] 



(4.28) 



where Yi,Y2,9 are independent random variables, Y,Yi,Y2 
have density /„, while 9 has density lia+i/2,a-i/2- The script 
has to be meant as an identity in distribution. To prove 



(14.28b it suffices to show that Y' ' 
alently. 



49[Yi+Y2\-^ or,equiv- 



X 



49X1X2 
X1+X2 



(4.29) 



where X , Xi , X2 are independent Gamma{a, 1 / (a — 1 ) ) random 
variables. The result follows if one is able to show that {Xi + 
X2){49XiX2y^ has Laplace transform (14271 1. Setting G 
Xi +X2 and B := Xi/{Xi +X2) one rewrites ( 14.29b as 



X ='U9GB{l-B). 



(4.30) 



It is a classical result of probability theory that, given Xi 
and X2 which are independent and Gamma distributed, the 
random variables G and B are stochastically independent 
and, moreover, that G has Gamma{2a,l/{a — 1)) distribu- 
tion, while B has Beta{a,a) distribution. For the proof of 
this property, we refer for instance to Chapter 10.4 in ifTsll . 
Moreover, using relation ( 13.111 ) one can reckon that 9G has 
Gamma (fl+l/2,l/(fl— 1)) distribution simply by computing 
its Laplace transform. Using now the fact that 9G and B are 
mutually independent, one can write the Laplace transform of 
49GB{1 - B) in the point £, as 



Jo n 



■ dx 



[l+4(fl-l)-i<^x(l-x) 



[l+4(fl-l)-i^.<l-x)]"+^ 



- dx. 
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At this stage, a simple change of variable shows that 



r(2fl) 



dz 



22"-ir(«)2 7o {l+^z/{a-\))Y+\ 

The last identity follows by the duplication formula r(2fl) — 
22"-' r(a)r (a + 1 /2) /r(l/2) (see 8.335.1 in |20]). Using re- 
lation ( 13.111 1 once again we get L((^) = (l+(^/(fl— 1))"" and 
( |4.29t is proved. 

Some remarks are in order Within the choice ( 14.141 ). the 
conservative in the mean trade ( 11.3b is such that the two agents 
maintain at least 1 /4 of the total wealth used to trade. Thus, 
trade ( 11.31 ) is in a sense less risky than the conservative trade 
( |1.2t . where one of the two agents can exit from the trade 
with almost no money. In addition, it follows that the number 
of moments of the explicit equilibrium state which are finite 
increase with a. On the other hand, when a increases, the area 
described by the distribution of the random fraction e on the 
interval [l,+oo) decreases, and the probability to use wealth 
of the society is also decreasing. Hence a fat Pareto tail is 
obtained through a big use of the common wealth. 



V. CONCLUSIONS 



where y/{x) = r'(x) /r{x) is the Digamma function ll20ll . Us- 
ing the dupUcation formula 2\(f{2x) = \ir(x) + \(f{x + 1/2) + 
21og(2), see 8.365.6 [20.] . one obtains 



G'(l) = -7 



1 r(2fl)r(fl-.s + i) 



4 T{2a - s)T{c. 



e(fl) 



where Q{a) := 21og(2) + \i/{a - 1/2) + \i/{a). 

Now 1/(1 /2) = —7— 21og2 and \jf{l) — -7 (7 being the 
Eulero-Mascheroni constant), see 8.366.1/2 |20], and then 

e(i)-o. 



The classical expansion formula (see 8.363.8 112011 ') 



k>0 



{x + kf 



allows to conclude that, for every a > 1 



In this paper, we introduced and discussed the equilibrium 
solution of a nonlinear kinetic equation of Boltzmann type, 
modelling redistribution of wealth in a simple market econ- 
omy in which trades are described by a standard gambling 
game. Due to the simplicity of the game trade, analytical so- 
lutions can be obtained in the case in which the post-trade 
wealths depend on the pre-trade ones through random vari- 
ables which are Beta distributed. Previously known analytical 
solutions are here shown to exit for particular values of the 
underlying parameters. Despite its simplicity, the model en- 
lightens the role of the interaction in producing Pareto tails. 



Appendix A: Appendix 



Let us prove that G'(l) < when G is defined as in ( 14.161 ). 
Starting from ( |4.16t , differentiation shows that 

G'(s) = _2^-2s'[!^fE(^^Zl±ll. 

r(2fl-s)r(a-t-i)' 

• {21og(2) + xi/{a-s + ^)- v/(2fl - i)} 



^'^""^ |i,(«-l/2 + fe)2 ^,{a + ky 



E 

k>0 ' 



>0. 



Hence Q{a) is strictly monotone in [l,+oo) and since Q{1) = 
it follows that Q{a) > for every a > 1. This shows that 
G'(l) <0. 
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